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In this paper we estimate the event rates of neutralino-nucleus scattering leading 
to electron emission by the atomic ionization. We find that the branching ratio 
for detecting electrons vis a vis the traditional neutralino nucleus elastic scattering 
sensitively depends on the threshold energy. In the case of a light target and a 
neutralino mass of 100 GeV we estimate it to be around 10 percent assuming a 
reasonable threshold energy of 250 eV. 
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1 Introduction 



The combined MAXIMA-1 [1], BOOMERANG [2], DASI [3] and 
COBE/DMR [4] Cosmic Microwave Background (CMB) obser- 
vations as well as the recent WMAP data [5] imply that the 
Universe is flat [6] and that most of the matter in the Universe 
is Dark, i.e. exotic. Crudely speaking one has: 

tt b = 0.05, ttcDM = 0.30, tt A = 0.65 

for the baryonic, dark matter and dark energy fractions respec- 
tively. Furthermore, since the non exotic component cannot ex- 
ceed 40% of the CDM [8], there is room for the exotic WIMP 
(Weakly Interacting Massive Particles). Many experiments are 
currently under way aiming at the direct detection of WIMP. 
In fact the DAMA experiment [9] has claimed the observation 
of such events, which with better statistics have subsequently 
been interpreted as a modulation signal [10]. These data, how- 
ever, are not consistent with other recent experiments, see e.g. 
EDELWEISS [11] and CDMS [12]. 

Supersymmetry naturally provides candidates for the dark mat- 
ter constituents [13], [14]- [17]. In the most favored scenario of su- 
persymmetry the lightest supersymmetric particle (LSP) can be 
simply described as a Majorana fermion, a linear combination of 
the neutral components of the gauginos and higgsinos [13], [14]- 
[24]. 

The event rates, however, are expected to be quite low and the 
nuclear recoil energies are extremely small. Thus one has to try 
to reduce the background to the lowest possible level and under- 
stand possible reaction induced backgrounds (for a brief discus- 
sion see the appendix at the end of the paper). Furthermore one 
has to search for characteristic signatures associated with this re- 
action. Examples are the modulation of the event rates with the 
motion of the Earth (modulation effect) and the correlation of 
the observed rates of directionally sensitive experiments with the 
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motion of the sun [25,26,27]. Transitions to low energy excited 
nuclear states have also been considered [28,29]. 

In the present paper we will explore another novel possibility, 
namely detecting dark matter constituents by observing the low 
energy electrons, which follow the ionization of the atom dur- 
ing the LSP-nucleus collision. This possibility may be realized 
with the technology of gaseous TPC detectors [30]. In fact the 
WIMP-nucleus scattering leads: (i) The nuclear recoil and (ii) to 
nuclear recoil with atomic excitation. So far most CDM searches 
have been made by the inclusive processes, (i) and (ii), employing 
solid detectors. We propose that the produced electrons in atomic 
excitation (ii) should be studied by exclusive measurements. 



2 The effect of bound electrons 



The differential cross section for the LSP nucleus scattering lead- 
ing to the emission of electrons in the case of non relativistic 
neutralinos takes the form: 

^4^I M | 2 (&S ( ^ )3 ^ 

E Pni [<Pnim(^r^(T x + € n , - T - ^- - (Px^"^ 



■ 2 2^, + ^-r-^-^- K - q 



where v, T x and p x are the oncoming LSP velocity, energy and 
momentum respectively, while q is the momentum transferred 
to the nucleus. M is the invariant amplitude, known from the 
standard neutralino nucleus cross section, T and k are the ki- 
netic energy and momentum of the outgoing electron and e n g is 
the binding energy of the initial electron. n ^ m (k) is the Fourier 
transform of the bound electron wave function, i.e its wave func- 
tion in momentum space. p n £ is the probability of finding the 
electron in the n, I orbit. 

In the expression above and in what follows our normalization 
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will consist of one electron per atom, to be compared with the 
cross section per nucleus of the standard experiments. 

The binding energies for hydrogenic wave functions are given by: 

_ m e {Za) 2 

€n ~ 2n 2 • [) 
One, of course, may have to introduce effective charges, espe- 
cially for heavy atoms. Some hydrogenic wave functions for light 
systems are going to be given below. 

After summing over the m-substates we find that: 

2(2^1) EI^(k)| 2 = ^fcW| 2 . (3) 

Due to the delta function the integration over the momentum q 
is trivial and yields 

J g dg6(T x + e nt -T- — — ) = vK^—^—, 

(4) 

with fi r the LSP-nucleus reduced mass and 



A = 

with 



tl ■ m X( 1 n m X (T-ent) 



Px Px 

In order to avoid any complications arising from questions re- 
garding the allowed SUSY parameter space, a subject extensively 
studied anyway [13]- [24], we will present our results normalized to 
the standard neutralino nucleus cross section. The thus obtained 
branching ratios are independent of all parameters of supersym- 
metry except the neutralino mass. The numerical results given 
here apply in the case of the coherent mode. If, however, we limit 
ourselves to the ratios of the relevant cross sections, we do not 
expect substantial changes in the case of the spin induced process. 
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With these ingredients we find that the ratio of the cross sec- 
tion with ionization divided by that of the standard neutralino- 
nucleus elastic scattering, nuclear recoil experiments (nrec) , takes 
the form: 



da(T) ir /n _ _ 2 



= -;J2Pni\(/>ni(2rn e T)l 

4 ni 



C"nrec 4 n g 

- (C+A) 2 

" " _ I / [//:,.( If I / 5. i I I 

-m e kdT, (7) 



Jo 1 2^[F(2//X)] 2 



where & = ^[1 + - 1)], 2^p x f = 2/^ is the mo- 

mentum g transferred to the nucleus and i^(<z) is the nuclear form 
factor. The outgoing electron energy lies in the range < T < 
-^-T x — e n £. The expression for the standard cross section a nrec 
will be given below (Eq. 9). 

Since the momentum of the outgoing electron is much smaller 
than the momentum of the oncoming neutralino, i.e. K « 1, the 
integration over £i can be trivially performed. Furthermore, if the 
effect of the nuclear form factor can be neglected, the integration 
over £ can be performed analytically. Thus we get: 

da(T) 1 ~ r- — - 2 m x (T-e n i) 
= oEM<Mv2m e T)| 7f\ ) 

@nrec ni f^r J- \ 
,m^{T — 6 n p) , , n = _ 



! _ ^S_jA )]me ^T dT , (8) 



Otherwise the angular integrations can only be done numerically. 
Furthermore integrating numerically the above expression over 
the electron spectrum we obtain the total cross section. One, of 
course, must convolute the above expression with the velocity 
distribution (see next section) to obtain both the differential rate 
as well the total rate as a function of the neutralino mass. 

Before proceeding further we will give an estimate of the expected 
standard elastic scattering recoil coherent events. The LSP- nu- 
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cleus cross-section takes the form: 



°nrec = (^) 2 A 2 <T p £ 2d^[F{2^)]\ (9) 

where o v is the LSP-proton cross section and fi r {p) ~ rn p is the 
reduced mass for the proton-LSP system. From this expression we 
see the coherent effect of all nucleons arising from the elementary 
scalar interaction. In the case of a light target like 20 Ne the effect 
of the form factor is not large. The event rate, corresponding to 
an LSP velocity \f{v 2 ), takes the form: 

R = § = P -^^4¥)(^?A^ (10) 

at m x Am p Hr\P) 

where p(0) = 0.3GeV/cm 3 is the LSP density in our vicinity and 
m is the detector mass. The upper limit of the current experi- 
ments is a p < 10 5 p6. Using this limit with m x = lOOGeV and 
\J{v 2 ) = 270 km I 's, we get for the A=20 system 

« 8 events/ (y. kg — target). 

This rate will be become « 10 events/ (y. kg — target), if the 
folding with the LSP velocity distribution is taken into account 
(see sec. 4). 



3 A brief discussion of the backgrounds 



We must begin by emphasizing that the standard dark matter ex- 
periments have thus far consisted of inclusive measurements em- 
ploying solid detectors, where the electrons and the recoil nuclei 
have not been separated. These measure continuum energy spec- 
tra of the recoiling nuclei with an energy resolution, which is not 
very good. Most experiments are employing neutrons with energy 
in the range of 5 — 10 MeV with a resolution of a few 0.1 MeV, 
which is much larger than the electron energy. Then the electron 
contribution is hidden. The present paper is concerned with an 
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exclusive measurement, whereby the electron signal is separated 
from the nuclear recoil signal. The electron rate is about 10% of 
the nuclear recoil rate. The actual nuclear recoil rate is reduced 
by about 30 — 50 % depending the detector ienergy cutoff. It 
thus appears that the electron rate could be as high as 1/3 of 
the nuclear recoil rate. This makes it a realistic proposal, if the 
BG problems are under control. The nice feature of an exclusive 
experiment is that one can measure simultaneously the electrons 
and the nuclear recoils, substantially reducing this way the BG 
rate. 

To circumvent the BG problems one must use high purity mate- 
rials and run the experiment underground. The main source of 
background is neutrons and this applies both in the standard ex- 
periments as well as the present novel proposal. It thus appears 
that background events of primary concern to the present work 
can arise from two main sources: 

(1) Reaction induced Auger electrons. These can arise from the 
decay of the inner electron holes in the daughter atom. The 
decay width has been estimated to be 1.2 x 10 _9 keV for 20 Ne. 
These electrons can have energies in the region of interest. In 
the case of 20 Ne they are expected at energies of 0.77, 0.80 
and 0.83 keV depending on which two electrons fill this hole 
and produce the Auger electron [30]. The reaction induced Is 
holes, however, occur with a small probability, see Fig. 1 One 
may also have background electrons, produced from impuri- 
ties which may lead to betas and photoelectrons producing 7 
rays. Anyway it has been estimated that such a background 
is under control [30]. 

(2) Radioactive Impurities (RI). RI 's give rise to (3, 7, a BG's. As 
it is well known, recent recoil detectors can separate the recoil 
signal from electron backgrounds of RI (3, 7's. Similarly the 
electron tracking detectors select electron signals. It is quite 
realistic, as shown both in the recent solar v detectors and 
the planned (3(3 experiments, to build detectors at the level of 
ppt of U-Th impurities (ppt= 10" 12 ). Note that most RI BG's 
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deposit (3, 7 energies at multi-sites. Then the BG rate for the 
electron detectors with impurities at the level of a few ppt is 
estimated at 30 per year per ton, by a stringent selection of 
DM single-site signals from the RI multi-site events. This is 
orders of magnitude less than the present limit (8/year/kg) 
of the WIMP's. It should also be noted that this limit is 
2-3 orders of magnitude less than the electron signal rates 
(1 ~ 10 per year per ton) and the electron BG rates of the 
planned electron detectors for real time solar p-p z/s and (3(3 
experiments. In any case it seems feasible to build detectors 
with BG rates well below the the signal rates with 10 _5 pb. 
Furthermore since the electron energy signal is always accom- 
panied by the recoil signal, the coincidence measurement of 
both the electron tracks and and the recoil signals will make 
this experiment almost free from BG events. 

(3) Electron backgrounds due to RI. These are under control [30] . 
Most electrons from RFs are (3 rays followed by 7 rays and 
Compton electrons followed by Compton 7 rays. These de- 
posit energy at multi-sites and, thus, they can be eliminated 
by selecting the single-site event from WIMP's. With a pu- 
rity level of the order of /iBq per Kg one gets a BG event 
rate of less than 0.1 per year per kg , after a stringent signal 
selection. This is far below the present WIMP signal rate (8 
per y per kg). 

(4) Neutron induced events. This source is also present in the 
standard nuclear recoil experiments and it should be dealt 
with in both. Such neutron events can be deciphered by 
studying the reaction products. We will, however, estimate 
the event rate due to neutrons using the distribution for muon 
produced neutrons given at 2000 m w.e. by Gaitskell [31]. 

• Neutron electron scattering due to the neutron magnetic 
dipole moment. 

This process is very special in the sense that the electron 
interacts directly with the neutron. Thus the magnetic mo- 
ment interaction does not suffer from suppression at high 
electron energies due to the bound electron form factor, see 
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Figs 4- 5 obtained as discussed in the appendix. 
• Electron ionization due to neutron-nucleus elastic scatter- 
ing. 

This process is characterized by kinematics similar to 
the reaction we are interested in. For the model discussed 
in the appendix involving a Yukawa interaction, see Fig. 6, 
we obtain the results shown in Figs 7-8. 
It appears that both of these backgrounds yield rates that 
could be troublesome. In practice the rate gets smaller by 
2-3 orders of magnitude at typical underground labs with 
4000 - 6000 m w.e. Active neutron shields for incoming and 
outgoing neutrons reduce further the BG rate. Accordingly 
the BG rate due to the neutrons can be negligibly small in 
realistic cases. These BG problems are further discussed in 
the Appendix. 



4 Some results 



We will now apply the above formalism for a typical case, namely 
m x = 100 GeV. We will employ a Maxwell-Boltzmann LSP veloc- 
ity distribution with respect to the galactic center [27], namely: 

/(v) = (^-<^> 

with vq = 220 km/s. On this we impose by hand an upper ve- 
locity bound (escape velocity), v esc = 2.84-uq. We then transform 
this distribution to the lab frame, v — > v + vq, where vq is the 
velocity of the sun around the center of the galaxy. We will not 
be concerned here with the motion of the the Earth around the 
sun, i.e. the modulation effect [27]. 

Folding both the numerator and the denominator of Eq. (7) with 
the LSP velocity distribution, after multiplying each with the 
LSP flux e^-f^v, we obtain the differential ratio with R e 

m x Am p ' R dT ' e 



the rate for ionization. If we then integrate the final expression 
with respect to the electron energy, we obtain the ratio By 
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doing the same in the case of Eq. (8), we obtain the analogous 
ratios and associated with the case without nuclear 

form factor. 

We will apply our results on a light target, e.g. 19 F, which is con- 
sidered in the standard nuclear recoil experiments. The obtained 
results are essentially identical to those for 20 Ne, which is a popu- 
lar gaseous TPC counter currently being considered for detection 
of low energy neutrinos produced in triton decay. Our approach 
can, of course, also be applied in the case of the heavier targets, 
employed in current dark matter searches, but the atomic physics 
involved is much more complicated 

The radial functions (f) n i(k) will be given in terms of a dimension- 
less variable b = by writing <j>ni(k) Rru(b), with R 2 nt {b) 
given by: 

2 256 
Rls = (1 + 46*)V (12) 

_ 512(l- 4 6 2 ) 2 _ 163846 2 

2s ~ (1 + U*)\ ' 2p ~3(l + 4&2)V 
The binding energies employed are given in the above order by 
[32]: 

e ls = -0.870 , e 2s = -0.048 , e 2p = -0.021. (14) 

In Fig.l we show the differential rate of our process, divided by 
the total nuclear recoil event rate, for each orbit (see Eq. 8) with- 
out the inclusion of the nuclear form factor on the left and with 
appropriate form factor [33] on the right. In Fig. 2 we show the 
same quantity taking all orbitals together. These results were ob- 
tained considering one electron per atom and weighing each orbit 
with the probabilities p n f. 

jto = (2/10, 2/10, 6/10). (15) 

The branching ratio for the total rate leading to our process, 
obtained as discussed above, are shown in Fig. 3, as a function of 
the low energy detector cut-off (threshold). From this plot we see 
that the introduction of the nuclear form factor, which we view 
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as quite reliable [33], does not substantially alter the branching 
ratio. This trend for the branching ratio may persist even for 
heavier nuclei, even though the individual rates may depend on 
the assumed form factor. 

From these plots we see that, even though the differential rate 
peaks at low energies, there remains substantial strength above 
the electron energy of 0.2 keV, which is the threshold energy used 
in MICRO MEG AS detectors, like the one recently [30] proposed. 



5 Conclusions 

Based on our results summarized in the previous section one can 
be optimistic about using the emitted electrons in the neutralino 
nucleus collisions for the direct detection of the LSP. This novel 
process may be exploited by the planned TPC low energy electron 
detectors. By achieving low energy thresholds of about 0.25 keV, 
the branching ratios are approximately 10 percent. They can be 
even larger, if one includes low energy cutoffs imposed by the 
detectors in the standard experiments, not included in the above 
estimate. 

As we have seen the background problems associated with the 
proposed mechanism are not worse than those entering the stan- 
dard experiments. In any case coincidence experiments with x- 
rays, produced following the de-excitation of the residual atom, 
may help reduce the background events to extremely low levels. 

Acknowledgments: This work was supported in part by the Euro- 
pean Union under the contracts RTN No HPRN-CT-2000-00148 
and MRTN-CT-2004-503369. One of the authors (JDV) is in- 
debted to Professors K. Zioutas for a careful reading of the manuscript 
and Y. Giomataris and N. Alexandropoulos for their useful com- 
ments regarding experimental issues related to this work. 
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Fig. 1. The differential rate, divided by the total rate associated with the nuclear 
recoils, as a function of the electron energy T (in keV). Each atomic orbit involved 
in the target 20 Ne is included separately. The results exhibited are obtained without 
nuclear form factor on the left and with nuclear form factor on the right. The full 
line, the short-dashed line and the long-dashed line correspond to the orbits Is , 2s 
and 2p respectively. 




0.4 T — > keV 




0.4 T — > keV 



Fig. 2. The same as in Fig. 1 by summing the contribution of all orbitals together. 




0.25 0.5 0.75 1 1.25 1.5 1.75 keV 




0.25 0.5 0.75 1 1.25 1.5 1.75 E±h _^ keV 



Fig. 3. The total rate for ionization divided by the standard nuclear recoil rate as 
a function of the threshold energy. No low energy cutoff imposed by the detector 
was included in the computation of the standard rate. The estimated rate for the 
proposed mechanism is obtained by multiplying this branching ratio with the stan- 
dard event rate. In the case of 20 Ne this about 10 events/ (kg-y) , obtained from the 
upper bound of the cross section as derived from the current experiments. 

Appendix: Some further discussion on the background. 



In this section we are going to briefly discuss the possible back- 
ground (BG) events. A more detailed discussion will be given 
elsewhere. RI impurities in all detector components should be as 
low as possible to avoid all kinds of BG electrons from (3 and 7 
rays from RI. In practice the Ur and Th impurities in the rock at 
underground laboratories can be reduced to be less than 1CT 16 in 
weight leading to a BG event rate of 0.03 per Kg per year, much 
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Fig. 4. The differential cross section for neutron electron scattering due to the 
neutron magnetic moment in the case 20 Ne. The results on the left correspond to 
an incident neutron energy of 2 keV plotted as a function of the outgoing electron 
energy for electrons bound in the Is, 2s and 2p orbitals. The results on the right 
correspond to an outgoing electron of 2 keVas a function of the incident neutron 
energy. 




Fig. 5. The differential event rate for electron ionization in neutron electron scatter- 
ing, which is due to the neutron magnetic moment, in the case 20 Ne as a function 
of the outgoing electron energy. Electrons bound in the Is, 2s and 2p orbitals were 
considered . The results presented were obtained using distribution of Gaitskell for 
muon produced neutrons. All neutrons were included, but the event rate reached 
saturation in about 500 MeV. 




Fig. 6. The strong interaction (Yukawa like) potential felt by the elastically scattered 
neutrons for the A = 20 system. 

lower than the DM event rate of the present concern. 

Neutrons in the 10 MeV -1 GeV region are of the order of 1 
per m 2 per year at underground laboratories of 4000m water 
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Fig. 7. The differential cross section for neutron electron scattering due to the strong 
interaction. The value a s = 1.0 was used in the numerical calculation. Otherwise 
the notation is the same as in Fig. 4. 




1 T^keV 



Fig. 8. The differential rate for electron ionization in elastic neutron nucleus scat- 
tering. Otherwise the notation is the same as in Fig. 5 

equivalent. Lower energy neutrons arising from nuclear fission of 
Ur-Th isotopes around the detector can be slowed down and be 
absorbed. 



This source is also present in the standard nuclear recoil exper- 
iments. Such neutron events can be deciphered by studying the 
reaction products. We will, however, estimate the event rate due 
to neutrons using the distribution for muon produced neutrons 
given by Gaitskell [31], approximated by:. 

f[e n ] = (i.7 2e -[ e «A 241MeV 0] + 8.27e-^ 55 - 5MeV ^)/MeV.y.m 2 . 

(A.l) 

We distinguish the following cases: 

• Neutron electron scattering due to the neutron magnetic dipole 
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moment. The leading term in this interaction is of the form: 

^ n r i2 n \s2/*\ 1 



Q(r) = V4^f^l(T e x a n } 2 ■ V^Y 2 (f)-. (A.2) 

Thus one must evaluate the orbital matrix element 

/ d 3 xQ(V2xy^J d'yfaC^jf), (A-3) 

where 4>b is the bound electron wave function and k is the mo- 
mentum of the outgoing electron. The evaluation of the above 
integral proceeds by a multipole expansion. For bound elec- 
trons in s-states the only multipole appearing is the quadrupole, 
while for p-states one can have both a dipole and octapole term 
with the dipole dominating. One can then compute the spin 
matrix elements and calculate the differential cross section as 
a function of the outgoing electron energy. One should use a 
proper neutron energy distribution function. For orientation 
purposes we calculated the differential cross section for a neu- 
tron energy of e n = 2 keV, which leads to a maximum electron 
energy e n + e n £. The obtained results are shown in Fig. 4. We 
see that the shape of the differential cross section is very dif- 
ferent from that of the interesting electrons. Furthermore such 
events can be rejected by the subsequent detection of the neu- 
tron. The event rate, obtained by using the Gaitskell neutron 
distribution, is given in Fig. 5. 

Neutron nucleus elastic scattering. This can lead to electron 
ejection in pretty much the same way as the LSP-Nucleus elas- 
tic scattering. Background events can be rejected by the detec- 
tion of the scattere neutron. 

Proceeding in a manner analogous to that described in sec. 
2 we find: 

da(k) = i^£| M |2 ^ 3 Pn A ( )3 I {AA) 
y J vE e l 1 (27r) 3 (27r) 3V ; 2(2^+1) v ; 

E Vnt [(^Wm(k)] 2 27r£(e n + e ne - T - e n - ^ Pn k — ^-), 

nim Zm A 

with the amplitude M being just the Fourier transform of the 
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potential generated by the A nucleons, which is felt by the 
neutrons: 

M = f dhe l{ ^-^- v V A {r) = V A (\p n - p n \). (A.5) 

Ignoring the energy taken away by the nucleus we finally find: 

I 

dcr(n,(e,n)) On, T) = — -^m 2 n m e V m e TdT (A.6) 

7TV Z 

^Pni\MV2^T)\ 2 F A (e, ^^), 
with F A (e n , x) given by: 



F A (e n , x) = J i VT _ tdt [V A {bV^h^ t)] , (A.7) 

where 6 is a length parameter characterizing the range of the 
potential. 

The essential input in this case is the neutron-nucleus elastic 
scattering potential V A , which can be determined phenomeno- 
logically. 

Proceeding further the elastic neutron nucleus cross-section 
can be obtained in an analogous fashion: 

d<7(„,„)(e n ) = f Q n ( e n ~ x)dx [V A (b(e n - x))] 2 . (A.8) 

In other words the evaluation of the interesting cross-section 
requires a complicated analysis of the elastic neutron nucleus 
cross-section. Thus for the purposes of the present work we 
have decided to obtain the potential by folding the Yukawa in- 
teraction with a range b = m^ 1 with a constant nuclear density 
. The thus obtained potential is shown in 6. The obtained elas- 
tic neutron nucleus rate is presented in Fig. 8. We see that this 
background contribution is peaked at very low electron energy. 
• Compound nucleus formation by neutron capture. Such back- 
ground events can be rejected by detecting the decay 7-rays of 
the compound nucleus. 

In this case, if one ignores the tiny energy taken away by 
the outgoing nucleus, for a neutron of energy e n the electron 
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energy is fixed, T = e n + e n £ + A — E x , where A is the available 
energy and E x is the energy of the of the compound nuclear 
state. For thermal neutrons the electron energy is essentially 
fixed by the energy of the compound nucleus T « e n £ + A — E x . 
Thus this source does not pose a serious problem. 

Anyway, if one wants, one can cast the relevant cross section 
in the form: 



2(27r) 2 \2e nn t; J 



p n e\/^rn e {e n + e ni + A - E x ) (A.9) 



\(j) ni {pm e {e n + e ne + A - E x ))\ 2 \M(p n , x)\\ 



where M(p n , x) is the amplitude associated with the standard 
compound nucleus formation. Noting that the captured neu- 
trons are thermal one finds that the event rate for this process 
can be related to the rate for neutron capture as follows: 



dN 1 re 

~17 ~ o/o ^ L de n^e E Pnt\l^rn e {e n + t n l + A - ^.10) 

at z^-KfJQ n ^ x 

\(J) ni {^2m e {e n + e ni + A- E x ))\' 



dt 



compound 



where the last term in the above equation gives the rate for 
the compound nucleus formation and e is the highest energy at 
which neutron capture can occur. We are not, however, going 
to elaborate further on this point. 

The background rates have so far been evaluated by using the 
neutron flux at 2000 m w.e. given by Gaitskell. The BG rates are 
of the same order of magnitude as the true DM rates or smaller 
and are well reduced by requiring anti-coincidence with the reac- 
tion products. Furthermore present underground laboratories are 
typically at around 3000 m w.e. or deeper. Thus the neutron flux 
and the background rates at these depths get smaller by orders 
of magnitude compared to those at 2000 m w.e. It thus appears 
that the BG rates from neutrons are not serious. 
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